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Abstract 

In this paper we study radial oscillations of the rotating strange 
stars in strong magnetic fields in the Density Dependent Quark Mass 
(DDQM) model. We see that increase of frequency i.e. difference in 
frequency of rotating and non-rotating stars is more for higher mag- 
netic fields. The change is small for low mass stars but it increases 
with the mass of the star. This change of frequency is significant for 
maximum mass whereas it is marginal for a 1AM@ star. 
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1 Introduction 



Since the original suggestion made by Comeron [1965] more than three decades 
ago that vibrations of neutron stars could excite motions that can have inter- 
esting physical applications, there have been several investigations of vibra- 
tional properties of neutron stars. It has been found that for typical neutron 
star models with mass about one solar mass and radius about lOkms these 
vibrations have periods in the range (3-5)ms and this period is relatively in- 
sensitive to the exact value of central density [Cutler, Lindblom and Splenter, 
1990]. 

The possibility that radial oscillations of neutron stars give rise to the 
oscillations observed with radio subpulses of pulsars was proposed by Bori- 
akoff [ 1976 ]. Inspite of doubts about this interpretation the possibility of 
detecting the radial oscillations have not been excluded [Van Horn, 1980]. 
More recently, periodicities have been observed in X-ray bursts [Sadeh et al, 
1982] which have raised considerable interest in radial oscillations [Glass and 
Lindblom, 1983; Muslinov and Tsygan, 1986; Marti et al, 1988]. Recently 
Vath and Chanmugan [1992] have determined the two lowest frequency ra- 
dial oscillation modes for several equations of state. They found that their 
results for neutron stars differ from those of Glass and Lindblom [1983] even 
if the central density is same in both the works. 

Rotation is a general property of all stellar bodies i.e. galaxies and in- 
dividual stars as evidenced by Doppler broadening of their spectral lines. 
When the core of a massive star collapses to produce neutron star which ul- 
timately may convert to a strange star [Witten, 1984; Farhi and Jaffe, 1984; 
Haensel, Zdunik and Schaeffer, 1986; Alcock, Farhi and Olinto, 1986 and 
Madsen and Haensel 1991] ; the conservation of angular momentum assures 
its enhanced rotation. Rotation has effects on structure and stability of rela- 
tivistic stars and on space-time. The rotation of local inertial frames affects 
the internal structure of rotating stars sometimes referred as dragging of the 
inertial frames or as the Lense-Thirring effect. The frame dragging frequency 
is largest at the centre of a rotating star and falls to zero at large distances 
where the space-time is flat. 

As rotation and magnetic field are intrinsic properties of all stellar objects, 
it is intresting to ask what will be the normal modes of a rotating magnetised 
star. Rotation can lead to coupling between oscillation modes. This can 
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influence the dissipation of radial vibrational energy of the star and also 
its rotational stability. The determination of normal mode frequencies of a 
rotating neutron star is non-trivial. The Newtonian formulation of this aspect 
of neutron stars has been studied only recently by Ipser and Lindblom [1989, 
1990]. Post-Newtonian frequencies for pulsations of rapidly rotating neutron 
stars have been reported by Cutler and Lindblom [1992]. The radial modes 
of rotating neutron stars in the Chandrasekhar-Friedman [1972] formalism 
have been studied by Datta et al [1998]. 

Since we are interested in studying the radial oscillation of magnetised 
rotating quark stars, we would like to mention that not much work has been 
done even for non-rotating magnetised quark stars. Recently Vath and Chan- 
mugan [1992], Gondek and Zdunik [1999], Benvenuto and Lugones [1998], 
Anand et al [2000a] have studied the radial oscillations of strange stars. Usu- 
ally the strange stars are studied in the framework of MIT Bag Model. An 
alternative description of quark matter in which the confinement is treated by 
assuming a baryon density dependence of the quark masses was introduced 
by Fowler, Raha and Weiner [1981] and Plumer, Raha and Weiner [1984]. 
Later on, the model was reformulated [Benvenuto and Lugones; 1995a, 1995b 
and 1998] to show that the properties of SQM in this model are quite similar 
to those predicted by MIT Bag Model. Anand et al (2000b) used the Den- 
sity Dependent Quark Mass (DDQM) model to study the bulk viscosity of 
SQM. DDQM model effectively includes the first order QCD coupling correc- 
tion and yet is much easier to work with. However in some respect DDQM 
model is quite different from the MIT Bag model. In the low pressure re- 
gion, the speed of sound is much larger than the relativistic c/ of the 
Bag model, and it agrees remarkably well with phenomenologically expected 
values from hadronic collisions over a large range of temperatures [Fowler et 
al(1981)]. Motivated by these arguments Anand et al [2000a] have studied 
the radial oscillation of magnetised quark stars in Density Dependent Quark 
Mass (DDQM) model. 

In this paper, we study the radial oscillations of magnetised rotating 
strange stars in DDQM model and compare the results with Anand et al 
(2000a) for non- rotating stars. The plan of the paper is as follows : in 
section 2, a brief review of the DDQM model is given. Section 3 deals with 
formalism and finally in section 4, the results and discussions are presented. 
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2 Density Dependent Quark Model in Mag- 
netic Field 



Following Benvenuto et al (1998) and Anand et al (2000a & b), we shall treat 
SQM as a free gas of u,d,s quarks and electrons where the mass of each quark 
is parametrised as 

m u = m d = C/3n B , m s = m so + C/3n B (1) 

Here m so , the strange quark current mass and C are arbitrary constants, 
restrained by the requirement that at T=0 

(a)E/n B < 930Mev for SQM 

and 



(b)E/n B > 940Mev (2) 

for two flavour quark matter. We have selected C = 75Mevfm~ 3 and m so = 
140Mev well in the stability window calculated by Benvenuto and Lugones 
(1998). 

We consider u, d, s, e system in the presence of a magnetic field B which 
is directed along the Z-axis. The energy of a spin 1/2 particle of mass rrii 
and charge qi in the presence of magnetic field is well known and is given by 



Si = \jm} + V \i + 2q t B(n + s + 1/2) (3) 

where s = +l/2(— 1/2) refer to spin up (down) states of the particle, and 
p Z) i is the momentum along the Z-axis. We have used /i/27r = c = 1. The 
thermodynamical potential of the system is given by 

n = En* = -E^E(2- <W / d Pz [in(i + e -^)\ - |vt 2 t 4 (4) 

where % = (u,d,s,e), gi = 6 for (w, d, s) and 2 for electrons. The last term 
in eqn. (4) is the contribution of gluons. The expressions for pressure pi, 
energy density and number density rii of various species are obtained by 
using the well known relations: 

Pi = {n B -^-) -fy (5) 
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E i = -p i + f i i n i -T&) (6) 

and 

711 = { d^ ] (7) 

T,n B 

At T=0, the above expressions reduce to 

g l q i B n -*\ 1 2 2m,C ^ + K fi 

Pi = -^-f Q (2-6 no )[- l MK fi -(m i +^B + ^)ln( - / —== )] 

(8) 

2^ ^ 2 2 3n B V«V + 2nqiB 

(9) 

and 

m = £ ( 2 - <W* (10) 

Z7r n=0 

where 

ify = ^2 _ m t 2 _ 2^5 (11) 

and 

rw* = mt[ (/ii2 ~^ 2) ] (12) 
2%-D 

As the system is in /^-equilibrium, we have 

lid = V* (13) 

and 

Hd = A*u + Me (14) 

(Assuming that neutrinos/antineutrinos produced stream out freely, fx v — 0). 
The charge neutrality condition gives 

2n u - n d - n s - 3n e = (15) 

and the baryon number density of the system is 

n B = (n u + n d + n s )/3 (16) 
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For a given baryon density and given magnetic field B the chemical poten- 
tials of various species are evaluated self-consistently by demanding charge 
neutrality and /3-equilibrium conditions. This enables us to make a profile 
of pressure (p = J2?i) and energy density (E = J2Ei) f° r various values of 
baryon densities and magnetic fields. (See Anand et al, 2000a). 



3 Radial Pulsations of a Slowly Rotating Mag- 
netised Quark Star 

We use the formalism given by Chandrasekhar and Friedman (1972) and 
Datta et al (1998) which is fully general relativistic to calculate the effect of 
rotation (to order Vt 2 , where Vt is the angular vlocity of rotation) on the eigen- 
frequencies of radial pulsations of stars. It gives an exact formula to calculate 
the frequency (a') of oscillations of a "slowly" rotating stellar configuration 
which depends only on the knowledge of the Lagrangian displacement associ- 
ated with the radial mode of oscillation of the non-rotating configuration and 
the uniform spherical deformation caused by rotation. By "slow" rotation 
we mean those rates of rotation that are small in comparision to the central 
mass shed limit. 

The Chandrasekhar-Friedman formula is of the form 



o' 2 I x = I 2 + I 3 + U (17) 



where 



h = J (fx^e 2X i 2 u \p + pc 2 ) (18) 
/ d 3 x^V 2 [ { ^ c2) (AN/N) 2 + 2rp(AN/N)±(ln Ul )Z 
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dr ' ->■■■■<■>■■' 

„2\„. „. ,-./ 



+ (P + ^ + 2(*V/V)*p + ^?(p + pc 2 ) 2 e 2 ^ 2 

5A 



+ 2(p + pc 2 )u u 1 Q , —l;] s (20) 



/ 



The subscript ,s, refers to Schwarzschild (i.e. non-rotating configuration) 
values and cu(r) is the angular velocity of the fluid element relative to the 
local inertial frame. 

In the above equations only the spherical deformation effects of rotation 
are manifest. Contributions due to quadrupole deformations are all taken to 
be zero. The metric used is similar to Hartle and Thorn (1968) i.e. 

ds z = - e 2 ^+^\cdt) 2 + r 2 sin 2 fl(# - udtf + r 2 d9 2 + e 2(x ° +x »\drf (22) 

where u(r) is the angular velocity of the commulative inertial frame dragging. 
The angular velocity of the local inertial frame is written as 

lu = fl - u (23) 

where Vt is the angular velocity as seen by a distant observer. Notice that 
for a non-rotating case the metric is given by 

ds 2 = -e~ 2v ° (cdtf + e 2As (drf + r 2 (d6 2 + sin 2 ^0 2 ) (24) 

In eqn (22) we have written 

v = v s + v Q (25) 

and 

A = X s + A n (26) 

where and \n are the corrections of O (Q 2 ) due to the rotations and v s 
and \ s are the Scharzschild values. 
In Hartle's notation 

m = K (27) 
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and 

Xn = r - 2Gm(r)/c 2 (2§) 
where h and m are functions of r. To o(f2 2 ), 

^ZT g = r 2 e (x+u) s [i + (A + v ) n ] sin 6 (29) 
r 2 e Us 

\f ZI g, = i sin 6 (30) 

1 r 2 u 2 

u = e -*(l -vn + -^e-^inH) (31) 
l c 

«i = rV sin = (r 2 cj/c)e^ s sin 2 # (32) 

V = (rio/c)e~ Ua sm6 (33) 
Various other quantities appearing in equations (19)-(22) are as follows: 

5\ = 5\ s + 5\ n (34) 
A-kG r(p + pc 2 ^ 

and 



^ = — ~ 2G^ ( 35 ) 



47rG rip + pc 2 ) r/ p + pc 2 2m r 2 oj 2 .,. 

«» = -— jhif^'H^^v^W^— e sm 9]( (36) 

AN/N = (AJV/JV), + (AJV/JV)„ (37) 



where 

<AW0. = £f-«A-f = -£ < 38 > 

where we have used £ = r~ 2 e Us ( 

(AN/N) n = [-^-(AN/N) S + (l/r)£]^e- 2 -sin 2 £ + ^ (39) 
p + per er ar 

^ = -[-^(AJV/JV)J-2J/r+^{-?{ (40) 
V p + pc 2 ar c<j 
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where 



dh a d 1 _ 2l; ra 2 

"dr^ 6 ( T } (41) 



4 d , 4 . da; N , ^ , dj 



1/r 4 — (r*j—) + = (42) 

' dr y dr' ' dr v 7 

dm Q 47rGV 2 de 1 j 2 r A dcu 2 1 qdj 2 _ 2 

— ^^ (e + ^ + i^ ( d^ "3? r d7" (43) 

^ m (l + ^g^) 47rG(6 + p)r 2 _L_ r 4 j 2 da) 2 

dr (r- 2 ^) 2 c 4 (r- P ° 12c 2 ( r - ^) 1 dr j 

+ — -( rW ) (44) 

The pressure and energy change from 
p — > p + Ap 
and 

e — »■ c 2 p + c 2 Ap. 
Hence, 

& | 2(v | + (P + Pc 2 ) 2 d(r p ) ( P + pc 2 ) 2 d Po 

dr C dr +Z[P + CP) Tp P ° (Tp) 2 dr P ° + (r P ) dr 

(46) 

For the sake of completeness we study the moment of inertia of a rotating 
strange star in strong magnetic field and the expression used is: 

/ = J/Q (47) 

where J, the angular momentum, is given by 

and 

n = q(r) + f§. 

The equation governing infinitesimal radial pulsations of a non-rotating 
star in general relativity was given by Chandrasekhar [1972] and has the 
following form 

F § + G fr +H ^^ < 48 > 
9 



where the coefficients F, G and H are given in detail by Anand et al [2000a] 
and Datta et al [1998]. The equations for hydrostatic equilibrium of non- 
rotating degenerate stars in general relativity are 



dp G(p + p/c 2 )(m + A-nr^p/c 2 ) 
dr = r 2 (l- ^) 

dm 



dv G(m + A7cr 3 p/c 2 
dr r 2 c 2 (l — 



(49) 



^ = 4vrr 2 p (50) 
dr 



(51) 



A = --Ml - — ) (52) 



1 2Gm 
'2 /n(1 -7c 

The boundary conditions to solve the pulsation eqns.(48) are 
£(r = 0)=0 
and 

5p(r = R) = -e| - rp~(r 2 e-»0\r=n = (53) 



4 Results and Discussion 

The calculation of the radial modes of vibration of a rotating star is straight- 
forward as long as the rotation is supposed to be "slow": il « il c = 
^GM/R 3 , but extremely tedious. For each star of the sequence correspond- 
ing to a given central density p c and magnetic field B, we have a profile of p, 
p, m, Tp, v s {r) etc as functions of r. The solution of the eigenvalue problem 
has provided us £,{r), £'( r ) for each mode. 

Our next task is to solve the second order differential equation for uJ c for 
a given cJ c (0) and obtain Q and J after the solution is obtained. The value of 
cJ c (0) is then adjusted to obtain each required value of Q. Next the equation 
for h Q (r) and p (r) are solved using the profile of a)(r) already obtained. Once 
this is done all the quantities involved in the integrals Ji, I2, I3 and I4 are 
known. The value of the frequency of radial oscillation a' for the rotating 
star can thus be found. 
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In fig. 1 we plot the moment of inertia I vs M/M & for both relativistic and 
non-relativistic strange stars for eB=0 and 10 5 Mev 2 . Here curves a and b 
correspond to non-relativistic case for eB=0 and 10 5 Mev 2 respectively while 
a', b' correspond to relativistic case for the same magnetic fields. It is seen 
that I increases with magnetic field for the same mass configuration for both 
non-relativistic and relativistic case. In fig.2 we have plotted uo/Vt vs r/R for 
M = M max i mum and for eB=0, 1 x 10 5 , 2 x 10 5 Mev 2 corresponding to a, b, 
c curves respectively. It is seen that the effect of rotation decreases as one 
moves from the centre to the surface of the star. However the magnetic field 
slows down the rotation considerably from the centre towards the surface 
of the star implying that the magnetised strange stars are less prone to the 
rotational instabilities. 

In figs. 3 to 5 we show the behaviour of a and a' (i.e. radial frequencies for 
non-rotating and rotating magnetised stars) for n=0, 1, 2 modes respectively 
for Q,=0.1il c . It is interesting to note that for non- rotating strange stars the 
radial oscillation frequencies decrease with increasing p c . For rotating stars 
the frequencies initially decrease but then start increasing. The increase 
is more pronounced for n=0 mode implying the increased stability of the 
magnetised strange stars. 

Figs. 6 to 8 show the variation of a and a' with M/M Q and depict a 
behaviour similar to the one in figs. 3 to 5. 

In figs.9-11, we show the plot of a' 2 - a 2 vs Q/Q c for eB=0, 1 x 10 5 , 
2 x 10 5 Mev 2 for n=0, 1, 2 modes respectively. It is noted that the increase in 
a' 2 — a 2 is monotonic with rotation and this effect increases with increasing 
magnetic field. 

To conclude we find that in the presence of a strong magnetic field the 
rotating strange stars are less prone to rotational instabilities implying they 
become more stable. 
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M/M e 



Figure 1: Moment of inertia, I, in units of M & Km 2 vs M/M Q for both 
relativistic and non-relativistic strange stars for eB = 0, 1 x 10 5 Mev 2 . Curves 
a and b correspond to the nonrelativistic case for eB = & 1 x 10 5 Mev 2 
respectively while a' and b' correspond to the relativistic case. 
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Figure 3: Radial oscillation frequencies a (non-rotating stars) and a' (rotat- 
ing stars) in units of 3 x 10 5 sec.^" 1 vs p c in units of 10 15 gm./cm. 3 . Curves 
a, b and c are for nonrotating stars while a', b' and d refer to rotating stars 
for magnetic fields 0, 1 x 10 5 and 2 x 10 5 Mev 2 respectively for n = and 

n = o.m c . 
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Figure 4: Same as figure 3 for n = 1. 
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Figure 9: (a' 2 — a 2 ) in units of 9 x 10 10 sec. 2 vs Q/Q c for n — 0. Curves a, 
b, c refer to maximum masses and corresponding magnetic fields: 2.17M , 
QMev 2 ; 3.OM , 1 x 10 5 Mev 2 and 3.66M , 2 x 10 5 Mev 2 respectively. 
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Figure 10: Same as figure 9 for n — 1. 
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Figure 11: Same as figure 9 for n = 2. 
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